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Abstract. 

We employ a time-dependent mean-field-hydrodynamic model to study the 
generation of bright solitons in a degenerate fermion-fermion mixture in a cigar- 
shaped geometry using variational and numerical methods. Due to a strong Pauli- 
blocking repulsion among identical spin-polarized fermions at short distances there 
cannot be bright solitons for repulsive interspecies interactions. Employing a linear 
stability analysis we demonstrate the formation of stable solitons due to modulational 
instability of a constant-amplitude solution of the model equations for a sufhciently 
attractive interspecies interaction. We perform a numerical stability analysis of these 
solitons and also demonstrate the formation of soliton trains by jumping the effective 
interspecies interaction from repulsive to attractive. These fermionic solitons can be 
formed and studied in laboratory with present technology. 
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1. Introduction 

After experimental observation and the study of bright sohtons in Bosc-Einstein 
condensates (BEC) [1-3], recent observations [4-7] and experimental [8-10] and 
theoretical [11-14] studies of a degenerate Fermi gas (DFG) by sympathetic cooling 
in the presence of a second boson or fermion component suggest the possibility of 
soliton formation [15] in a degenerate fermion-fermion mixture (DFFM). Apart from 
the observation of a DFG in the degenerate boson-fermion mixtures (DBFM) ^'^Li [6], 
^^Na-^Li [7] and ^'''Rb-^°K [8,9], there have been studies of the following spin-polarized 
DFFM ^oK-4°K [4] and ^Li-^Li [5]. 

The dimensionless one-dimensional nonlinear Schrodinger (NLS) equation in the 
attractive (self-focussing) case [16] 

d(f) 1 d'^d) , ,,9 , „ /X 

sustains the following bright soliton [16]: 

t)^a sech[a(y - yt)Yvy-i{v^-a^)t/2+i^ (2) 

The parameter a represents the amplitude as well as pulse width, v the velocity, and cr is 
a phase constant. These bright solitons are possible only in one dimension. In two and 
three dimensions they are allowed in the presence of a transverse trap [1, 15], or in the 
presence of an oscillating nonlincarity [17]. Apart from bright solitons in the attractive 
case, gap solitons are possible in the repulsive (self-defocussing) case in the presence of 
a periodic potential [18]. 

Bright solitons in a BEC are formed due to a nonlinear atomic attraction [1,2]. As 
the interaction in a pure DFG at short distances is repulsive due to strong Pauh blocking, 
there cannot be bright sohtons in a DFG. However, bright sohtons can be formed [19,20] 
in a DBFM in the presence of a sufficiently strong boson-fermion attraction which can 
overcome the Pauli repulsion among identical fermions. Bright solitons can also be 
formed in a binary mixture of repulsive bosons supported by interspecies attraction [21]. 

We demonstrate the formation of stable fermionic bright solitons in a DFFM for 
a sufficiently attractive interspecies interaction. In a DFFM, the coupled system can 
lower its energy by forming high density regions, the bright solitons, when the attraction 
between the two types of fermions is large enough to overcome the Pauli repulsion. 
We use a coupled time-dependent mean-field-hydrodynamic model for a DFFM and 
consider the formation of axially-free localized bright solitons in a quasi-one-dimensional 
cigar-shaped geometry using numerical and variational solutions. The present model 
is inspired by the success of a similar model used recently in the investigation of 
collapse [14,22] and bright [20] and dark [23] solitons in a DBFM, and black sohtons [24] 
and mixing- demixing [25] in a DFFM. 

We study the condition of modulational instability of a constant-amplitude solution 
of the present model under a plane-wave perturbation and demonstrate the possibility 
of the formation of bright sohtons by a hnear stabihty analysis. We find that for 
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a sufficiently strong interspecies interaction, under the plane-wave perturbation the 
constant-amplitude solution becomes unstable and localized solitonic solutions may 
appear. We present a numerical stability analysis of these bright solitons by introducing 
different small perturbations when the solitons undergo stable and sustained breathing 
oscillation. 

We also consider and study the formation of a soliton train in a DFFM by a 
large sudden jump in the interspecies fermion-fermion scattering length realized by 
manipulating a background magnetic field near a Feshbach resonance, experimentally 
observed in two hyperfine states of the following spin-polarized DFFMs: ^Li-^Li and 
40j^_40j^ [26]. This procedure effectively and rapidly turns a repulsive DFFM to a 
highly attractive one and generates bright solitons. An experimental realization of 
bright solitons in a DFFM could be tried in two hyperfine states in samples such as 
®Li-®Li or ^"K-^'^K using a Feshbach resonance. However, there is already experimental 
evidence [27] and theoretical conjecture [28] that the maximum attractive force, that 
can be created by a Feshbach resonance in fermionic atoms in two hyperfine states is 
limited by quantum mechanical constraints of unitarity. Although there is a limit to 
the creation of attraction in a two-component spin-polarized DFFM in two hyperfine 
states, there is no such limit in a multi-component spin-polarized DFFM [29] or a DFFM 
composed of atoms of distinct mass [30] . So, if bright solitons and soliton trains cannot 
be efficiently created in a two-component DFFM in two hyperfine states, a better and 
more efficient DFFM for the creation of solitons could be the mixture of fermionic 
atoms of distinct mass, where one can avoid the problem of a possible suppression of 
interspecies attraction. One such system, among many others, could be the mixture 
of spin-polarized ^Li-'^^'K mixture: both ®Li [5, 31] and '^'^K [4] have been trapped and 
studied in laboratory. The formation of bright solitons in a DFFM by turning the 
effective interspecies interaction among spin-polarized fermions from repulsion to strong 
attraction seems within the reach of present experimental possibilities. 

Here we shall be interested in the formation of bright solitons in a spin-polarized 
DFFM of different fermionic atoms in the presence of strong interspices attraction and 
we shall ignore the possibility of the formation of a Bardcen-Cooper-Schreiffer (BCS) 
condensate [32] through Cooper pairing. A BCS condensate is usually formed in the 
presence of a weak attraction among identical fermions with fermion pairs in the singlet 
(spin parallel) state. The formation of a BCS condensate should not be favored in a 
spin-polarized (spin antiparallel) two-component DFFM with strong attraction among 
the different types of atoms, as a strong interspecies attraction is not the domain of BCS 
condensation. By choosing a strong interspecies attraction our study stays in the BEC 
region of the BCS-Bose crossover problem [33] strongly favoring molecule formation. 
The solitons once formed will decay via molecule formation of two different types of 
fermionic atoms [34] as in the case of bosonic solitons. Nevertheless, such molecule 
formation is a slow process and can be accounted for by a three-body recombination 
term as introduced in a previous study of collapse in a DFFM [35]. The molecule 
formation will eventually destroy the solitons in the DFFM, albeit, at a slow rate. The 
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same is true in the formation of soliton and soliton train in an attractive BEC, where the 
essential features of the dynamics have been well explained within the mean-field Gross- 
Pitaevskii model by neglecting molecule formation [36]. A subsequent study including 
the effect of molecule formation [37] has not essentially changed the conclusions of [36]. 
Hence this pioneering study on the formation of soliton and soliton trains in a DFFM of 
different atoms using mean-field hydrodynamic equations with the neglect of molecule 
formation is expected to explain the essential features of the dynamics. However, it 
would be of interest to include the effect of molecule formation on the DFFM solitons 
in a future study. 

In section 2 we present our mean-field-hydrodynamic model and its reduction to 
a quasi-one-dimensional form in a cigar-shaped geometry. In section 3 we show that 
bright solitons can appear in this model through modulational instability of a constant 
amplitude solution. In section 4 we perform a variational analysis of the mean-field 
equations. Numerical results for isolated bright solitons are presented in section 5 and 
are compared with variational results. Then we study the generation of a train of solitons 
by modulational instability. Finally, we present a summary in section 6. 

2. Nonlinear Model 

We use a simplified mean-field-hydrodynamic Lagrangian for a DFG used successfully to 
study a DBFM [14,20,23]. To develop a set of time-dependent mean-field-hydrodynamic 
equations for the interacting DFFM, we use the following Lagrangian density [14, 20] 



where j = 1,2 represents the two components, ipj the probability amplitude, nj = 
the probability density, * denotes complex conjugate, rrij the mass, Aj = 
^^(67r^)^/^/(2mj), the interspecies coupling gi2 = 27ih'^ai2/mR with mji = mim2/{mi + 
7712) the reduced mass, and ai2 the interspecies fermion-fermion scattering length. The 
number of fermionic atoms Nj is given by / drnj{r) = Nj. The trap potential with 
axial symmetry is taken as Vj{r) = ^Smju'^^p^ + u'^z^) where u and uu are the angular 
frequencies in the radial (p) and axial (z) directions with u the anisotropy. The u 
limit corresponds to a cigar-shaped geometry and allows a reduction of the three- 
dimensional equations to a quasi-one-dimensional form appropriate for freely moving 
solitons. The interaction between identical intra-species fermions in spin-polarized state 
is highly suppressed due to Pauli blocking terms 3Ajn^j^^/5 and has been neglected in 
(IH]). The kinetic energy terms fi^|Vr'?/'jP /{Grrij) in ([3]) contribute little to this problem 
compared to the dominating Pauli-blocking terms. However, its inclusion leads to an 
analytic solution for the probability density everywhere [20]. 




2 



+ 



(3) 
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With the Lagrangian density ([3]), the following Euler-Lagrange equations: 

d dC ^ d dC dC , , 



dt ^ dxk dtpi* 

^ dt k=l " ^ dxk 

with Xk, k = 1,2,3 being the three space components, become 



= 0, (5) 



where j ^ k = 1,2. This is a time-dependent version of a similar time- independent 
hydrodynamic model for fermions [13]. For large nj, both lead to [13, 14] the Thomas- 
Fermi result [3] Uj = [(/ij — Vj)/Aj]^^'^ with jjj the chemical potential. They yield 
identical results for time-independent stationary states. However, the present time- 
dependent model can be used in the study of nonequilibrium dynamics, as in the study 
of soliton trains. 

For z/ = 0, ([5]) can be reduced to an effective one- dimensional form by considering 
solutions of the type il)j{r,t) = \/N j(j)j{z,t)il)f\p) where 

I ,(0)/ _ ( mujp^ \ 

1^, (p)| =^exp^-^J, (6) 

with m = Srrij corresponds to the ground state wave function in the radial trap alone in 
the absence of nonlinear interactions. Here to have an algebraic simplification we have 
taken the masses of two types of fermions to be equal (mi = m2). These wave functions 
satisfy 

- ^ VJV^]°) + Imco'p'^f = h^i^f , (7) 

with normalization 27r {p)^^ pdp = 1. Now the dynamics is carried by t) and 

the radial dependence is frozen in the ground state ip^^\p). In the quasi-one-dimensional 
cigar-shaped geometry the linear fermionic probability densities are given by \(l)j{z,t)\'^ . 

Averaging over the radial mode il)f\p), i.e., multiplying ([5]) by il)f'**{p) and 
integrating over p, we obtain the following one-dimensional dynamical equations [25]: 



(8) 



where 



Njmu 



In dH]) the normalization is given by \(j)j{z,t)\'^dz = 1. In these equations we have 
set the anisotropy parameter z/ = to remove the axial trap and thus to generate 
axially-free quasi-one-dimensional solitons. 
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To reduce three-dimensional equations (IHl) to a dimensionless form, following [20], 
we consider variables r = tuj/2, y = z/l, (fi = \A(f)i, with I 
becomes 

92 



>yh/ (cum), while 



.d_ 



(9) 



where Njj 



9(67riVj)2/75, and Njk 



corresponding to attraction, and / 



00 
00 



12|ai2|A^i!c//. Here we employ a negative 012 
iPj{y,T)\'^dy = 1. In ([9]) a sufficiently strong 
attractive fermion-fermion coupling Njk\(fk\'^{j 7^ k) can overcome the Pauli repulsion 
Njjl^pjl^^^ and form bright solitons. 

For the usual Gross-Pitaevskii with the cubic nonlinearity, the reduction of the 
full three-dimensional equation to its one-dimensional counterpart was performed, in 
different forms, by many authors [38]. The nonlinearity with 4/3 power in was 
obtained starting from a similar nonlinearity in a three-dimensional DFG. A strict 
one-dimensional (two-dimensional) DFG will give rise to a quintic (cubic) nonlinearity. 
Whatever be the nonlinearity in the intraspecies fermions, bright fermionic solitons will 
be generated, provided that the interspecies attraction is strong enough to overcome the 
intraspecies repulsion due to Pauli blocking. However, in this paper we shall consider 
only the Pauli-blo eking nonlinearity with the 4/3 power as in (Q. 

The two coupled equations ([9]) could be simplified for Ni = N2 = N, while 
ipi = ip2 = ip, and these equations reduce to the following single equation 



0, 



(10) 



- dr dy"^ 

where [3 = Nu = N22 and 7 = A''i2 = A^2i- Equation ffTUl) maintains the essential 
features of (Q, e.g., a quadratic nonlinear attraction and a Pauli blocking repulsion. 
The one-dimensional equation (fTOj) is quite similar in structure to the NLS equation ([1]) 
apart for the Pauli-blocking repulsive term f3. For a small (3 and large 7 the solitons of 
([T]) survive in (fTOl) . However, they disappear in the opposite limit of large (3 and small 

7- 



3. Modulational Instability 

3.1. Symmetric Case (Ni = N2) 

We find that ffTUj) allows a constant-amplitude solution which exhibits modulational 
instability leading to a modulation of the solution. We perform a stability analysis of this 
solution and study the possibility of generation of solitons by modulational instability 
in the symmetric case. We consider the constant-amplitude solution [16] 

^0 = Aoexp{t5) = Aoexp[ti-fAlT - PA'J^)] (11) 

of ffTOj) . where Aq is the constant amplitude and 6 a phase. The time evolution of solution 
( fTTj) maintains the constant amplitude Aq but acquires an amplitude-dependent phase. 
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Now we study if this solution is stable against small perturbations by performing a 
linear stability analysis. 

We consider a small perturbation of the constant-amplitude solution ( ITTil given by: 

ip=iAo + A)exp{i6), (12) 

where A = A{y, r) is the small perturbation. Substituting the perturbed solution (|T^ 
in (ITOl) . and for small perturbations retaining only the linear terms in A we get 

dA d'^A 2 



(13) 



We consider the complex plane-wave perturbation 

A{y, t) = Ai cos{Kt - ny) + iA2 sin(Kr - Qy) (14) 

in ffT^ . where Ai and A2 are the amplitudes of the real and imaginary parts, respectively, 
and ii' is a frequency parameter and Q a wave number. Then separating the real and 
imaginary parts we get 

-AiK = A2n\ (15) 

-A2K =Ain^~2^AlAi + ^(3Al^'Au (16) 

and eliminating Ai and A2 we obtain the dispersion relation 

(17) 



12 '^/?/l4/3^ 



- (27^^ - -/3Al 



K = ±Q 

The constant-amplitude solution (JTTi) is stable if perturbations at any wave number Q 
do not grow with time. This is true as long as frequency K is real. From f|T7|) we 
find that K remains real for any Q provided that 2'yAl < 4(3 A^^^/ 3 or 7^0^^ < 2(3/3. 
However, K can become imaginary for 7^9^^ > 2(3/3 and the plane- wave perturbations 
can grow exponentially with time r. This is the domain of modulational instability of 
a constant-intensity solution [16]. The perturbation then grows exponentially with the 
intensity given by the growth rate or the modulational instability gain g{VL) defined by 

g{n) ^ 2^iK) = 2\n\ [2^Al - pA^ - fi^] (18) 

where 3" denotes imaginary part. The presence of modulational instability is closely 
connected with the appearance of a bright soliton [16]. Localized bright solitons are 
possible only when the constant-amplitude solution is unstable. 

3.2. Asymmetric Case (Ni 7^ N2) 

Now we consider the possibility of modulational instability of a similar constant-intensity 
solution in coupled equations We consider the constant-amplitude solutions 

iPjQ = Ajoexp{i6j) = Ajoexp[iT{NjkA^jQ - NjjAfo^)], 
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of ([9]), where Ajo is the amphtude and 6j a phase for component j. The constant- 

amphtude solution develops an amplitude dependent phase on time evolution. We 
consider a small perturbation Ajexp{i6j) to these solutions via 

(fj = {Ajo + Aj) exp{i5j), (19) 

where Aj = Aj{y,T). Substituting this perturbed solution in IQ, and for small 
perturbations retaining only the linear terms in Aj we get 



^ + (20) 



+ NjkAkoAjoiAk + Al) = 0, k. 

We consider the complex plane-wave perturbation 

Aj{y, t) = Aji cos{Kt - ny) + iAj2 sm{KT - Qy) (21) 

in fl20p for j = 1,2, where Aji and Aj2 are the amplitudes for the real and imaginary 
parts, respectively, and K and Q are frequency and wave numbers. Then separating the 
real and imaginary parts we get 

- AuK = Au^^ (22) 



-AuK = An^^ - 2N^2AioA2oA2i + ^NuA%^Au, 



(23) 



for j = 1, and 

- A21K = A22n^ (24) 



-A22K = A2in^ - 2N2lAwA2oAu + -A^224o'^2i, 



for j = 2. Eliminating A12 between fl22p and (1231) we obtain 

Au[K^ - + iNnA%^] = 2A2iNi2AwA2on\ 



(25) 



(26) 



and eliminating A22 between (HM and ( l25i) we obtain 

A2i[K^ - n^n^ + 4:N22A%''] = 2AiA^2i^io^2ofi'- 

(27) 

Eliminating and A21 from fl^Uj) and flTTj) . finally, we obtain the following dispersion 



relation 



= ±n[(n' + '^NuA%' + In22A%') ± (28) 



±1^ [NnAf - N22Ati'y + AN,2N2iAlAiy^' 



1/2 
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For stability of the plane-wave perturbation, K has to be real. For any this 
happens for 

or for N^2N2lA%^Al^'' < 4iVniV22/9. However, for N^2N2lA%''A%^ > AN^^N^i/Q [39], 
K can become imaginary and the plane-wave perturbation can grow exponentially with 
time. This is the domain of modulational instability of a constant-intensity solution 
signalling the possibility of coupled fermionic bright soliton to appear. In the symmetric 
case Ni = N2 and = A20 = Aq, consequently, Nn = N22 = P and N12 = N21 = 7 

2 /3 

and we recover the condition of modulational instability ^ 2/?/3 derived in section 

3.1. 



4. Variational Analysis 

4-1- Symmetric Case (Ni = N2) 

Next we present a variational analysis of (fTOl) based on the Gaussian trial wave function 
[40] 

,2 



(Pv{y,T) = Aexp 



y 



(29) 



where A is the amplitude, R is the width, b the chirp, and c the phase. The Lagrangian 
density for (fTUj) is the one-term version of ([S]), e.g.. 



C = -h 
2 



~ dt 



dip] dip 



+ 



dy 



which is evaluated with this trial function and the effective Lagrangian L = C{ipy)dy 
becomes 



^^AW^.6v^ 7^2 



(30) 



+ 



1 



b^R' 



2 i?2 

where the overhead dot denotes time derivative. The variational Lagrange equations 
d dL _dL 
dt dq dq ' 

where q stands for c. A, R and b can then be written as 



(31) 



A^R = = constant. 



(32) 
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2 ,2 .V^A^/' 4c 

--26^ + 2v/2,--4-/5^--. 



36 



2 

2/?6. 



66^ + V2' 



7 



i?2 



12^3 
5 



4c 



(33) 

(34) 
(35) 



The constant in (!32l) is fixed by the normahzation condition. Ehminating c from ([33 
and flM|) we obtain 



■21 4^3 



(36) 



The use of §2) 
R: 



351) and ( l36l) leads to the following differential equation for the width 



(PR 



d 

dR 



2 



2 

R 



1 8/5^3 \ 
2 1 12/?V3 



(37) 



(38) 



i?2/3 5vri/3v/5J' 

The quantity in the square bracket is the effective potential of the equation of motion. 
Small oscillation around a stable configuration is possible when there is a minimum in 
this potential. The variational result for width R follows by setting the right hand side 
of (1371) to zero corresponding to a minimum in this effective potential, from which the 
variational profile for the soliton can be obtained [40] . 



4-2. Asymmetric Case (Ni ^ N2) 

The above variational analysis can be extended to the asymmetric case. However, the 
algebra becomes quite involved if we take a general variational trial wave function with 
chirp and phase parameters. As we are interested mostly in the density profiles, we 
consider the following normalized Gaussian trial wave function for fermion component 
jof © 

' ,2 



1 



R,{t)^ 



exp 



y 



J 



1,2. 



(39) 



Using essentially the Lagrangian density 
Lagrangian 



in this case we obtain the following effective 



A^,- / 6^3 



+ 



with j ^ k = 1,2. The variational Lagrange equations ( 1311) for Ri and R2 now become 
4 4iV,, R, , 8v^ N,, _ ^ ^^^^ 



R^ 



+ 



[Rl + Riy/^ 5v^7rV3^5/3 



Equations (1401) can be solved for variational widths Rj and consequently the variational 



profile of the wave functions obtained from (]39l) . When Ni = N2, in (1401) N. 
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Figure 1. The solitons \(t)j{z)\ of © vs. z for iVi = 40, N2 = 60, ai2 = -300 
nm, wliile Nn « 149, iVi2 = 216, iV2i = 144, and N22 ~ 195. The corresponding 
variational solutions (j)jy{z) calculated from the coupled equations ()40p are also shown. 



Njk = Nkj = 7 and Ri = R2 = R] and in this case fj^Ul) yields the same variational 
widths as from the result obtained in the symmetric case in section 4.1 given by ( 1371) : 



5. Numerical Results 

We solve ([9]) for bright solitons numerically using a time-iteration method based on the 
Crank-Nicholson discretization scheme [41]. We discretize the coupled partial differential 
equations ([9]) using time step 0.0002 and space step 0.015 in the domain —8<y<8. 
The second derivative in y is discretized by a three-point finite-difference rule and the 
first derivative in r by a two-point finite-difference rule. We perform a time evolution 
of ([9]) introducing an harmonic oscillator potential y'^ in it and setting the nonlinear 
terms to zero, and starting with the eigenfunction of the harmonic oscillator problem: 
^i{y, t) = TT^^/^ exp(— ?/^/2) exp(— ir). The extra harmonic oscillator potential, set equal 
to zero in the end, only aids in starting the time evolution with an exact analytic form. 
With this initial solution we perform time evolution of ([9]). During the time evolution the 
nonlinear terms are switched on slowly and the harmonic oscillator potential is switched 
off slowly and the time evolution continued to obtain the final converged solutions. In 
addition to solving the coupled equations (Q, we also solved the single equation ( ITOl) in 
the symmetric case: A^i = N2. For given values of A^i and A^2; solitons can be obtained 
for 1 012! above a certain value. In the coupled-channel case solitons are easily obtained 
for a smaller |ai2| when A^^i and N2 are not very different from each other. 




(41) 
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In our numerical study we take / = 1 jjm and consider a DFFM consisting of 
two electronic states of atoms. This corresponds to a radial trap of frequency 
u = h/{l'^m) 27r X 83 Hz. Consequently, the unit of time is2/uj ^ A ms. For another 
fermionic atom the u value gets changed accordingly for / = 1 fim. 

5.1. Single Soliton 

From a solution of coupled equations ([9]) we find that these equations permit solitonic 
solutions provided that |ai2| is larger than a certain threshold value consistent with the 
analysis of section 3. First we solve ([9]) for Ni = 40, = 60, and au = —300 nm. In this 
case no solitons are allowed for ai2 = —290 nm. The solitonic solution suddenly appears 
as 1 012! increases past 290 nm. The solitons in this case are shown in figure 1, where we 
also plot the variational solutions of coupled equations fHOj) . In this case the nonlinearity 
parameters are A^n ^ 149, A^i2 = 216, N21 = 144, and N22 ~ 195. Substituting 
these values of the nonlinearities in ( l40l) and solving we obtain the variational widths 
Ri ~ 0.043 and R2 ~ 0.050. From fl29l) and we then obtain the variational soliton 
profiles \(f)vi{z)\ ^ 3.62 exp(-270z2) and \(f)v2{z)\ ^ 3.36 exp(-200z2). From figure 
1 we find that the variational results agree well with the numerical solutions of the 
coupled equations. In this case we also found the variational result in the symmetric 
case (A^i = N2 = 50) using (HTl) . which lies close to the above two variational solutions. 

To test the robustness of these solitons we inflicted different perturbations on them 
and studied the resultant dynamics numerically. First, after the formation of the solitons 
of figure 1 with A^i = 40 and N2 = 60 we suddenly changed the particle numbers to 
A^^i = 60 and N2 = 40 at time t = 100 ms. This corresponds to a sudden change of 
nonlinearities from A^n ^ 149, A^i2 = 216, A^2i = 144, and N22 ^ 195 to A^n ^ 195, 
A^i2 = 144, A^2i = 216, and N22 ~ 149. The resultant dynamics is shown in figures 2a 
and 2b. Due to the sudden change in nonlinearities, the fermionic bright solitons are 
set into stable non-periodic small-amplitude breathing oscillation. This demonstrates 
the robustness of the solitons. 

After the formation of the solitons of figure 1 with A^^i = 40 and A'^2 = 60 we 
suddenly changed the interspecies scattering length ai2 from —300 nm to —330 nm at 
time t = 100 ms. This can be realized by manipulating a background magnetic field 
near a fermion-fermion Feshbach resonance [26] This corresponds to a sudden change 
of nonlinearities from A^n ^ 149, A^i2 = 216, A^2i = 144, and A^22 ~ 195 to A^n ^ 149, 
A'^12 ~ 238, A"2i ~ 158, and A^22 ~ 195. Due to the sudden change in nonlinearities, 
the fermionic bright solitons are set into stable non-periodic small-amplitude breathing 
oscillation. Instead of plotting the soliton profile in this case we plot the root mean 
square (rms) size of the solitons in figure 3 which demonstrates the stable nonperiodic 
breathing oscillation. We also (i) gave a small displacement between the centers of these 
solitons and (ii) suddenly changed 0i — 1.1 x 0i and 02 1-1 x 02- In both cases after 
oscillation and dissipation the solitons continue stable propagation which shows their 
robust nature. 
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Figure 2. The propagation of fermionic solitons (a) |(^i(z,t)| and (b) \(l)2{z,t)\ of 
figure 1 vs. z and t for A^i = 40 and N2 = 60. At t — 100 ms (marked by arrows) 
the bright solitons are set into small breathing oscillation by suddenly changing the 
fermion numbers to A^i = 60 and N2 =40. 

5.2. Soliton Train 

During the time evolution of if the nonhnearities are changed by a small amount 
or changed slowly, usually one gets a single stable soliton when the final nonhnearities 
are appropriate. However, if the interspecies attraction is increased suddenly by a 
large amount by jumping ai2 from a positive (repulsive) value to a large negative 
(attractive) value, a soliton train is formed as in the experiment with BEC [1] because 
of modulational instability [42]. 

To illustrate the formation of soliton train in a fermion-fermion mixture in numerical 
simulation, we consider the solution of ([9]) for Ni = N2 = 40 and ai2 = —300 
nm with an added axial harmonic trap y"^ corresponding to nonhnearities Njj ^ 149 
and Njk ~ 144, j ^ k = 1,2. After the formation of the solitons in the axial trap 
we suddenly jump the scattering length to —400 nm, corresponding to off-diagonal 
nonhnearities Njk = 192, and also switch off the harmonic trap at time t = 0. Although 
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Figure 3. The rms sizes of tlie solitons and |(/i2(2,t)| of figure 1 vs. t. At 

t = 100 ms the solitons of figure 1 are set into small breathing oscillation by suddenly 
changing ai2 from —300 nm to —330 nm. 



the initial value of scattering length (ai2 = —300 nm) in this case corresponds to an 
interatomic attraction, because of strong Pauli-blocking repulsion the effective fermion- 
fermion interaction is repulsive in this case. However, the final value of scattering length 
(ai2 = —400 nm) corresponds to a stronger interatomic attraction which can overcome 
the Pauli repulsion so that an effective fermion-fermion attraction emerges in this case 
which might allow the formation of soliton(s). In our numerical simulation we find that 
this is indeed the case. Upon a jump of the scattering length to ai2 = —400 nm, we find 
that after some initial noise and dissipation the time evolution of ([9]) generates a single 
stable bright soliton as shown in figure 4a. 

However, more solitons in the form of a soliton train can be formed for a larger 
jump in the scattering length. In figure 4b, from the same initial state in figure 4a we 
consider a larger jump of to —570 nm corresponding to off-diagonal nonlinearities 
Njk ~ 274. The final value of scattering length {au = —570 nm) in this case corresponds 
to a stronger interatomic attraction than considered in figure 4a which might allow the 
formation of soliton trains. This is verified in numerical simulation. Upon a jump of 
the scattering length to = —570 nm, from ai2 = —300 nm, we find that after some 
initial noise and dissipation the time evolution of ([H]) generates three slowly receding 
bright solitons as shown in figure 4b. More solitons can be generated when the jump in 
the scattering length ai2 or off-diagonal nonlinearities is larger. 

In figure 4c we show the generation of five receding solitons of each component 
upon a sudden jump of the scattering length to —750 nm from the initial state of 
figure 4a. This corresponds to a jump of the off-diagonal nonlinearities to Njk ~ 360. 
The formation of soliton trains from a stable initial state is due to modulational 
instability [16]. The sudden jump in the off-diagonal nonlinearities could be effected by 
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Figure 4. Soliton trains of one, tliree and five solitons formed for Ni — N2 — 40 
upon removing the liarmonic trap and jumping the nonlinearities at t = 



from Ni 



7V.„ 



U9,Njk ~ U4:,k ^ j = 1,2 to (a) A^^ 



149, iV, 



jk 



192, (b) 



149, Njk ~ 274, and to (c) Njj « 149, Njk ~ 360, respectively, corresponding to 
a jump in scattering length ai2 from —300 nm to —400 nm, —570 nm, and —750 nm. 
For this purpose we solved the coupled equations Q. 
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a jump in the interspecies scattering length au obtained by manipulating a background 
magnetic field near a Feshbach resonance [26]. 

6. Summary 

We used a coupled mean-field-hydro dynamic model for a DFFM to study the formation 
of bright solitons and soliton trains in a quasi-one-dimensional geometry by numerical 
and variational methods. We find that an attractive interspecies interaction can 
overcome the Pauli repulsion and form fermionic bright solitons in a DFFM. This is 
similar to the formation of bright solitons in a coupled boson-boson [21] and boson- 
fermion [19, 20] mixtures supported by interspecies interaction. We show by a linear 
stability analysis that when the interspecies attraction is larger than a threshold value, 
bright solitons can be formed due to modulational instability of a constant-amplitude 
solution of the nonlinear equations describing the DFFM. 

The stability of these solitons is demonstrated numerically through their sustained 
breathing oscillation initiated by a sudden small perturbation. We also illustrate 
the creation of soliton trains upon a sudden large jump in interspecies attraction by 
manipulating a background magnetic field near a Feshbach resonance [26] resulting in 
a sudden jump in the off-diagonal nonlinearities. This jump transforms an effectively 
repulsive DFFM into an effectively attractive one, responsible for the formation of a 
soliton train due to modulational instability. Bright solitons and soliton trains have 
been created experimentally in attractive BECs in the presence of a radial trap only 
without any axial trap [1] in a similar fashion by transforming a repulsive BEC into an 
attractive one. In view of this, fermionic bright solitons and trains could be created in 
laboratory in a DFFM in a quasi-one-dimensional configuration. 

Here we used a set of mean- field equations for the DFFM. A proper treatment 
of a DFG or DFFM should be done using a fully antisymmetrized many-body Slater 
determinant wave function [11, 19, 43] as in the case of scattering involving many 
electrons [44]. However, in view of the success of a fermionic mean-field model in studies 
of collapse [14], bright [20] and dark [23] soliton in a DBFM and of mixing and demixing 
in a DFFM [25], we do not believe that the present study on bright solitons in a DFFM 
to be so peculiar as to have no general validity. 
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